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ABSTRACT 

We present a new homogeneous sample of 32 spectroscopic binary orbits in the 
young (~ 150 Myr) main-sequence open cluster M35. The distribution of orbital 
eccentricity vs. orbital period (e — log(P)) displays a distinct transition from 
eccentric to circular orbits at an orbital period of ~ 10 days. The transition is due 
to tidal circularization of the closest binaries. The population of binary orbits in 
M35 provide a significantly improved constraint on the rate of tidal circularization 
at an age of 150 Myr. We propose a new and more robust diagnostic of the degree 
of tidal circularization in a binary population based on a functional fit to the 
e — log(P) distribution. We call this new measure the tidal circularization period. 
The tidal circularization period of a binary population represents the orbital 
period at which a binary orbit with the most frequent initial orbital eccentricity 
circularizes (defined as e = 0.01) at the age of the population. We determine the 
tidal circularization period for M35 as well as for 7 additional binary populations 
spanning ages from the pre main-sequence (~ 3 Myr) to late main-sequence (~ 
10 Gyr), and use Monte Carlo error analysis to determine the uncertainties on 
the derived circularization periods. We conclude that current theories of tidal 
circularization cannot account for the distribution of tidal circularization periods 
with population age. 
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1. INTRODUCTION 

In a coeval population of binary stars the closest binaries (separations < 0.5 AU) tend 
to have circular orbits, while wider binaries have orbits with a distribution of non-zero eccen- 
tricities. This trend holds true for binaries with early- type (North & Zahn 2003; Pan et al. 
1998; Matthews & Mathieu 1992; Giuricin et al. 1984) as well as late-type (Mathieu et al. 
2004; Latham et al. 2002; Melo et al. 2001; Mathieu et al. 1992) main-sequence components, 
for binaries with evolved stellar components (MermiUiod & Mayor 1992), and in the orbital 
eccentricity distribution of large planets aroTind solar-type stars (http://cxoplancts.org). 
The transition from eccentric to circular orbits is explained by dissipative tidal interactions 
caused by the reaction of the binary components to each other's gravitational field (tidal cir- 
cularization: Hut 1981; Zahn 1977; Darwin 1879). The theory of tidal interactions predicts 
that the timescale of circularization depends strongly on stellar separation. Consequently, 
for a coeval binary population, the transition from eccentric to circular orbits should occur 
at a well-defined binary period. However, despite the coeval nature of a star cluster, its 
population of short period binaries may be small in number and consist of systems with a 
distribution of initial eccentricities, masses, and stellar angular momenta. Accordingly, the 
size and heterogeneity of the binary population will add complexity and uncertainty to the 
observational determination of the tidal circularization period. 

Nonetheless, the transition between eccentric and circular binary orbits provides an 
important constraint on the rate of tidal circularization as a function of orbital period, in- 
tegrated over the lifetime of a binary population. Furthermore, the distribution of tidal 
circularization periods with population age enables us to study the evolution of tidal circu- 
larization with time. 

The success of theoretical models describing the efficiency and the evolution of tidal 
circularization can be measured by their ability to "predict" the tidal circularization periods 
for binary populations of different ages. The differences between present-day models lie in 
the mechanism by which energy and angular momentum is transported between the binary 
components and their orbits. Currently, three theoretical models exist for the mechanism for 
tidal dissipation in late-type main-sequence stars: 1) The equilibrium tide theory describes 
retardation of the hydrostatic tidal bulge (equilibrium tide) due to the coupling of the tidal 
flow to the motion of turbulent eddies in the stellar convective envelope. This dissipative 
couphng is assumed to be responsible for a phase shift between the tidal bulge and the 
orbital motion of the binary stars. Because of the phase shifts a tidal torque is established 
between the two stars (Zahn 1977; Hut 1981); 2) The theory of dynamical tides describes 
the excitation and damping of gravity (g) waves in the radiative zones of stars due to the 
tidal forcing by the companion star (Savonije & Papaloizou 1983; Zahn 1977, 1975, 1970). 
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This dissipation mechanism has been successfully applied to explain circularization of early- 
type main-sequence stars with radiative envelopes (Claret & Cunha 1997; Giuricin et al. 
1984; Zahn 1977), and recently has also been applied to the radiative cores of late-type 
stars (Savonije & Witte 2002; Witte & Savonije 2002; Goodman & Dickson 1998; Terquem 
et al. 1998). Inclusion of resonance locking between the tidally forced modes and stellar 
eigenmodes provides increased efficiency of the tidal coupling (Savonije & Witte 2002; Witte 
& Savonije 2002, 2001, 1999a,b); 3) A pure hydrodynamical mechanismhas been proposed by 
Tassoul (2000, 1988) and Tassoul & Tassoul (1996). Tassoul suggests that hydrodynamical 
flows (large-scale meridional currents) induced by lack of symmetry about the rotation axis 
in the tidally perturbed star are responsible for the tidal torques on the component stars. 
This mechanism has been controversial on theoretical grounds (Rieutord 1992; Rieutord & 
Zahn 1997; Tassoul & Tassoul 1997). 

Despite the developments in the theory of tidal circularization, the models still cannot 
account for the observed periods of circular orbits in solar-mass binary populations (e.g. 
Witte & Savonije 2002; Claret & Cunha 1997; Mathieu et al. 1992). To help along further 
progress of theoretical modeling there is a need for high quality observational data on large 
and coeval samples of binary stars, allowing for accurate determination of the transition 
from eccentric to circular orbits. 

This paper presents a homogeneous {Mpnm ~ 1 ^o) sample of 32 spectroscopic binary 
orbits in the open cluster M35 (NGC 2168; ctsooo = 6^ 9™, ^2000 = 24° 20'; / = 186.59°, 
6 = 2.19°, distance ~ 850 pc). With an age similar to the Pleiades but approximately three 
times richer, M35 is a benchmark for stellar astrophysics at ~ 150 Myr (Deliyannis et al. 
2004, in preparation). The age makes M35 an important testing ground for binary stars that 
have recently ended their pre main-sequence (PMS) phase where they were larger and fully 
convective. Thus M35 sets the initial state for main-sequence tidal circularization. 

Section 2 describes the observations, instrument, and telescope used in the spectro- 
scopic survey of M35. In Section 3 we present the distribution of orbital eccentricity as a 
function of orbital period in M35 (the e — log(P) diagram). In Section 4 we motivate and 
present a new robust method for determining the period at which the most frequent binary 
orbits circularize. We will call this the tidal circularization period. In Section 5 we discuss 
and compare our new diagnostic to the previous measure of tidal circularization, the tidal 
circularization cutoff period. Wc compare the precision and accuracy performance of both 
diagnostics and determine measurement uncertainties on the tidal circularization period. In 
Section 6 we determine the tidal circularization period for the binary population in M35, and 
in Section 7 we determine the circularization period for a series of published populations of 
late-type unevolved main-sequence binaries. In Section 8 we carry out a comparison between 
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current theoretical predictions of the efficiency and evolution of tidal circularization and the 
circularization periods derived in Sections 6 and 7. Section 9 summarizes and presents our 
conclusions. 



2. OBSERVATIONS 

M35 has been included in the WIYN Open Cluster Study (WOCS) program since 1997. 
The radial-velocity survey of the cluster will be described in detail in Meibom et al. (in 
preparation). We summarize the relevant points here. The initial selection of ~ 2000 target 
stars was based on photometric membership in the cluster and proper-motion membership 
studies to V < 15 by McNamara & Sekiguchi (1986) and Cudworth (1971). Bjorkman 
& Mathieu (unpublished) completed astrometry and photometry to K ~ 17. Our radial- 
velocity study includes stars from y ~ 13 to V ~ 16.5, corresponding to a range in stellar 
mass from ~ 1.4 Mq {{B - V)o ~ 0.37) to ~ 0.7 Mq {{B - V)o ~ 1.1), with solar mass stars 
at V ~ 15 ((B-V) ~ 0.58). Masses are derived using the Yale (Y2) stellar evolution models 
(Yi et al. 2003). The cluster reddening {E(^b-v) — 0.20) was adopted from Deliyannis et al. 
2004 (in preparation). 

All spectroscopic data were obtained using the WIYN^ 3.5m telescope at Kitt Peak, Ari- 
zona, USA. The telescope is equipped with a Multi- Object Spectrograph (MOS) consisting 
of a fiber optic positioner (Hydra) feeding a bench mounted spectrograph. The Hydra posi- 
tioner is capable of placing ~ 95 fibers in a 1-degree diameter field with a precision of 0.2". 
In the field of M35 approximately 82-85 fibers are positioned on stars while the remaining 
fibers are used for measurements of the sky background. We use the 3" diameter fibers opti- 
mized for blue transmission, and the spectrograph is configured with an echelle grating and 
an all-transmission optics camera providing high throughput at a resolution of ~ 20,000. All 
observations were done at central wavelengths of 5130A or 6385 A with a wavelength range of 
~ 200A, providing rich arrays of narrow absorption lines. Radial velocities with a precision 
of ~ 0.5 kms~^ (Meibom et al. 2001) are derived from the spectra via cross-correlation with 
a high S/N sky spectrum. 

Telescope time granted from NOAO^ and Wisconsin allowed for 3-4 observing runs per 



■^The WIYN Observatory is a joint facility of the University of Wisconsin-Madison, Indiana University, 
Yale University, and the National Optical Astronomy Observatories. 

^NOAO is the national center for ground-based nighttime astronomy in the United States and is operated 
by the Association of Universities for Research in Astronomy (AURA), Inc. under cooperative agreement 
with the National Science Foundation. 
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year for M35, with each run typically including multiple observations on several sequential 
nights. Once identified, velocity variables are observed at a frequency appropriate to the 
timescale of their variation. At present the radial-velocity survey has resulted in a sample of 
81 spectroscopic binaries: 32 members, 37 candidate members and 12 non-members. Of the 
32 member binaries, 25 are single-lined (SBl) and 7 are double-lined (SB2) spectroscopic 
binaries. The orbital periods span from 2.25 days to 1115 days among the member bina- 
ries, corresponding to separations from ~ 0.04 to 2.5 AU assuming a 1 Mq primary and a 
0.5 Mq secondary component. The 32 member binaries were identified in the 2MASS All Sky 
Survey based on equatorial coordinates. 2MASS ID's, astrometry, photometry, and orbital 
parameters are listed in Table 1. 



Table 1. Astrometric, photometric and orbital data for binary members of M35. 
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3. THE PERIOD-ECCENTRICITY DISTRIBUTION IN M35 

The new compilation of binary orbits in M35 greatly improves the constraint on tidal 
circularization at the early stage of the main-sequence phase. Prior to this publication, the 
small sample of binary orbits in the Pleiades provided only rough limits on the transition 
between eccentric and circular orbits at this age (Mermilliod et al. 1992). Figure 1 shows 
the distribution of orbital eccentricity (e) as a function of log orbital period (P) for the 32 
spectroscopic binary members of M35. The shortest period binaries have orbital eccentricities 
consistent with circular orbits, while the orbits for longer period binaries show a distribution 
of non-zero eccentricities. The presumably primordial distribution of eccentricities for the 
longest period binaries is consistent with the Gaussian distribution observed in late-type 
main-sequence binary populations in other open clusters (Duquennoy et al. 1992), in the 
solar neighborhood (Duquennoy & Mayor 1991), and from the Galactic halo (Latham et al. 
2002). The lack of high eccentricity orbits (e > 0.6) among the long period binaries in M35 
may be due in part to observational biases. 

All binaries whose periods and eccentricities are plotted in Figure 1 are cluster members 
based on their radial velocities and location in the color- magnitude diagram (CMD). Fifteen 
have confirmation of membership from proper- motion measurements. A few binary systems 
deserve special attention and Figure 2 shows on a hnear scale the distribution of orbital 
eccentricity vs. period for binaries with periods less than 20 days. The last 4 digits of the 
2MASS ID are used to label key binaries in Figure 2 and for reference in the text and in Table 
1. Nine binary orbits have eccentricities less than 3 x cTe above e = 0.0, where cXe denotes the 
error on the orbital eccentricity. We consider these binary orbits to be circular. Binary 4037 
[Porb = 16.49 days, e = 0.041 ±0.014) has the longest period among the binaries with circular 
orbits. Three circular binaries have orbital periods close to 10 days: 6089 {Porb = 10.08 days, 
e = 0.008 ± 0.005), 0447 {Porb = 10.28 days, e = 0.009 ± 0.019), and 6200 (Porb = 10.33 
days, e = 0.016 ± 0.009). 4 binaries have eccentric orbits with periods shortward of the 
longest period circular orbit: 0422 (P = 8.17 days, e = 0.538 ± 0.013), 3081 (P = 12.28 
days, e = 0.550 ± 0.003), 7000 (P = 12.57 days, e = 0.095 ± 0.006), and 6415 (P = 14.52 
days, e = 0.246 ±0.044). 

Thus regardless of the homogeneity in age, metaUicity and primary mass of the M35 
binary population there is a ~ 8 day overlap in period between eccentric and circular binary 
orbits. 
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Fig. 1. — The distribution of orbital eccentricity (e) as a function of log orbital period (P) 
for 32 solar- type spectroscopic binary members in M35. 
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Fig. 2. — The eccentricity distribution of M35 binary orbits with periods shorter than 20 
days. Numbers next to selected data refer to binary ID numbers in Table 1. 
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4. DETERMINING THE CIRCULARIZATION PERIOD 

The transition from eccentric to circular orbits observed in M35 is also found in obser- 
vational studies of late-type main- sequence binary populations in well known open clusters 
such as the Pleiades (Mermilliod et al. 1992), the Hyades/Praesepe (Duquennoy et al. 1992; 
MermiUiod & Mayor 1999), M67 (Latham et al. 1992), and NGC188 (Mathieu et al. 2004), 
and in populations of PMS (Melo et al. 2001, and references therein), field (Duquennoy & 
Mayor 1991), and halo (Latham et al. 2002) binaries. This characteristic fingerprint of tidal 
circularization is an important observational constraint on the theory of tidal circularization. 
Consequently, determination of the period at which binary orbits becomes circular in a given 
binary population is critical to constraining theoretical models. However, determining this 
tidal circularization period is not straightforward in any of the above populations. 

The longest period circular orbit has been the preferred measure over the past decade 
(e.g. Duquennoy et al. 1992) and has been referred to as the tidal circularization cutoff' period 
(hereinafter: cutoff period). We argue here that the cutoff period is not the optimal measure 
of the transition between eccentric and circular orbits. 

Populations of short period binaries vary in size and consist of systems with a dis- 
tribution of initial eccentricities. This is observed in the Gaussian distribution of orbital 
eccentricity at longer orbital periods where presumably orbital evolution has been minimal 
(Halbwachs et al. 2003; Duquennoy et al. 1992; Duquennoy & Mayor 1991). The Gaussian 
distribution of initial orbital eccentricities combined with the effects of tidal circularization 
will lead to an overlap in period-space between initially low eccentricity orbits that have 
been circularized and initially high eccentricity orbits not yet circularized. This overlap is 
indeed observed in all published populations of coeval binary orbits with the exception of the 
small sample in the Pleiades cluster (see Section 7, Figure 8). The width of this overlap was 
estimated by Mathieu & Mazeh (1988) based on the observed distribution of initial orbital 
eccentricities, and indeed, Duquennoy et al. (1992) showed that the existence of eccentric 
orbits in the Hyades/Praesepe populations with periods shorter than the longest period cir- 
cular orbit can be explained due to incomplete tidal circularization of binaries with initially 
high orbital eccentricity. We similarly find through simulations of tidal circularization using 
the equilibrium tide theory by Zahn (1977) that eccentric orbits with periods shorter than 
the longest period circular orbit can be a result of high initial eccentricity (e > 0.5). 

For short-period binaries with low but non-zero orbital eccentricity (0.0 < e < 0.1), 
Mazeh (1990) suggests that the eccentricity is induced by the presence of a third companion 
star due to dynamical interactions different from tidal effects. However, most eccentric orbits 
with periods shorter than the longest period circular orbit have eccentricities between 0.1 
and 0.6, and among the 4 binaries in M35 with moderate to significant eccentricities and 
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periods shortward of the longest period circular orbit we see no indications in the stellar 
spectra of a third component. 

Monte Carlo simulations of tidal circularization using the equilibrium tide theory by 
Zahn (1977) show that the longest period circular binary orbit defining the cutoff period is 
likely to originate from the low eccentricity tail of the initial Gaussian eccentricity distri- 
bution. Accordingly, the cutoff period does not measure tidal circularization of the most 
frequent binary orbit. We show in Figure 3 the e — log(P) diagram of a tidally circularized 
population of artificially generated binaries (grey diamonds). To focus on the process of 
tidal circularization, the lower limit on the initial orbital eccentricity was set to 0.05 (dotted 
horizontal line). The threshold of 0.05 is 2-3 times the typical error on orbital eccentricities 
determined from observations. The initial and final locations of the longest period circular 
orbit are marked with an open and a filled circle, respectively. 

Figure 3 also demonstrate that the probability of measuring the "true" cutoff period 
of a parent population of binaries is highly sensitive to the size of the observed binary 
population. Consequently, the value of the observed cutoff period will vary with the size of 
the population. 

Furthermore, long period (P ~ 2 x 10^ — 10^ days) circular (e < 0.05) or marginally 
circular (e < 0.1) orbits of unevolved binaries are observed in the majority of the published 
binary populations (see Section 7, Figure 8) , indicating the possibility that binary stars can 
also form with very low eccentricity. While not included in Figure 3, the cutoff period is 
vulnerable to such primordial circular orbits. 

Finally, the determination of the cutoff period makes no use of the information provided 
by eccentric orbits in the e — log(P) diagram. 

Consequently we suggest that the cutoff period is not a robust measure of the state of 
tidal circularization in a coeval binary population. 

So motivated, we propose a new method for determining the period of circularization of 
primordial eccentric binary orbits. We refer to this period as the tidal circularization period. 
Our proposed method will make use of the information provided by all binary orbits in a 
population by fitting to the observed period-eccentricity distribution a function of the form: 




(1) 



This mathematical function is motivated by the observed period-eccentricity distribu- 
tions and by numerical modeling using the theory of Zahn (1977). The function is not 
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Fig. 3. — The e — log(P) diagram for a tidally circularized population of artificially generated 
binaries (grey diamonds). The lower limit on the initial orbital eccentricity is shown as a 
dotted fine ai e — 0.05. The initial and final locations of the longest period circular orbit 
are marked with an open and a filled circle, respectively. 



- 14- 

physically derived but constructed to mimic the tidal circularization isochrone of the most 
frequently occurring eccentric binary orbits. Figure 4 shows the function in the e — log(P) 
diagram. The transition of the function from eccentric to circular orbits is managed by the 
qI3{p -p) term and the 7 coefficient. 7 controls the abruptness of the break from e = at 
the period P' and f3 controls the overall steepness/slope of the transition. Our choices of f3 
and 7 are discussed in Section 5. The value of a (a = 0.35) is set to ensure that e(P > P') 
approaches a value of 0.35, the mean eccentricity of all observed binary orbits with periods 
longer than 50 days in the Pleiades, M35, Hyades/Praesepe, M67, and NGC188. For periods 
shortward of P' the function is set to 0.0 (e(P < P') = 0.0). We will hereafter refer to this 
function as the circularization function. 

The value of P' determines the location of the onset of the rise in the circularization 
function. We determine the location of the circularization function in period by minimiz- 
ing the total absolute deviation (S) between the observed (cj) and the calculated (e(Pj)) 
eccentricities, 

TV 

d = J2\ei- e{Pi) r . (2) 
1=1 

N denotes the number of binary orbits. The value of the exponent rj controls the influence 
on P' by the binary orbits with maximum deviation from the circularization function. We 
found that rj > 1 causes high sensitivity to the short period eccentric binaries, often leading 
to values of P' smaller than the period of the shortest period eccentric orbit, rj — 1 leads to 
values of P' between the periods of the shortest period eccentric orbit and the longest period 
circular orbit. We use rj — 1 when fitting the circularization function throughout this paper. 

Figure 4 shows the circularization function fitted to an artificially generated sample of 
binary orbits. Tidal circularization of these binary orbits has been simulated by numerical 
integration of differential equations derived using equations (4.3) and (4.4) in Zahn (1977) 
(Note erratum by Zahn (1978)) and the assumptions by Duquennoy et al. (1992): 

de —e 
dt ^ APW^ 

dP_ _ -3e^ 

The constant A controls the effectiveness of tidal dissipation and thus the rate of tidal 
circularization. We set A so that numerical simulation of tidal circularization using equations 



(3) 



(4) 
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(3) and (4) produce transitions from eccentric to circular orbits at periods similar to those 
observed. 

The circularization function mimics the transition from eccentric to circular orbits of 
the "typical" binary (0.3 < ej„j < 0.4; highlighted in Figure 4). The tidal circularization 
isochrone for the typical binary is reproduced by the circularization function at e = 0.01. 
Thus we define the "tidal circularization period" (hereinafter circularization period or CP) as 
the period for which the best fit circularization function has a value of 0.01 (e(CP) = 0.01). 
The choice of e = 0.01 as the threshold for circularization is also similar to the threshold cho- 
sen in theoretical simulations (see Section 8), and thus facilitates direct comparison between 
observations and theory. 

Our estimate of the initial eccentricity of a binary orbit that circularizes with a period 
equal to the circularization period at the age of a binary population, is tied to the tidal theory 
of Zahn (1977). Deviations from the typical (most frequent) initial eccentricity depends on 
the value of the circularization period (the age of the binary population). By fitting the 
circularization function to rich distributions of tidally circularized binary orbits (as shown in 
Figure 4) , we find that over the range of circularization periods of interest to this study, ~ 5 
days to 15 days, the range of initial eccentricities of the binaries represented by the resulting 
circularization periods is 0.2-0.5. This range of initial orbital eccentricities is centered on 
the mean of the observed Gaussian eccentricity distribution (e = 0.35) and we conclude that 
according to the Zahn (1977) theory the circularization period represents the orbital period 
at which a binary orbit with the most frequent initial orbital eccentricity circularizes. The 
range of initial eccentricities from 0.2 to 0.5 also includes the initial eccentricity used in 
theoretical predictions of tidal circularization (see Section 8) facilitating direct comparison 
between these predictions and the observed circularization periods. 



5. EVALUATION OF THE PERFORMANCE OF THE TIDAL 

CIRCULARIZATION PERIOD 

In the following we test the performance of the circularization function and compare its 
performance to the cutoff period. Using a Monte Carlo approach, we generate 10,000 binary 
populations each with 20 binaries with periods between 10 and 50 days, corresponding to 
the number of known binary orbits in the M35 sample with orbital periods shorter than 
50 days. Tidal circularization was simulated by numerical integration of the Zahn (1977) 
differential equations (eqs. [3] and [4] above). Initial orbital periods and eccentricities for 
each population were determined by random selection from a log-normal period distribution 
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Fig. 4. — The circularization function (black solid line) fitted to an artificially generated 
sample of binary orbits (diamonds). Tidal circularization of these binary orbits has been 
simulated by numerical integration of the differential equations (3) and (4). Selection of 
the initial orbital periods and the initial orbital eccentricities is described in the main text. 
Binaries with the "typical" (most frequent) initial eccentricities between 0.3 and 0.4 are 
highhghted (dark diamonds). The circularization period (CP) is defined as the period at 
which the circularization function has a value of 0.01 (e(CP) = 0.01). e = 0.01 is marked 
by a horizontal dotted line and CP is marked by the vertical black line on the period axis. 
The insert is a close-up view of the transition-region. 
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and a normal (Gaussian) eccentricity distribution. ^ We adopted the period distribution 
of Duquennoy & Mayor (1991) derived from solar-type binaries in the solar neighborhood, 
while the Gaussian eccentricity distribution (e = 0.35,0" = 0.21) was determined by fitting 
the distributions of all orbital eccentricities from the Pleiades, M35, Hyades/Praesepe, M67, 
and NGC188 for orbits with periods longer than 50 days. 

We ran the Monte Carlo experiment twice. The first time we excluded initially circular 
orbits (cj < 0.05) to test the performance of the diagnostics in the absence of initially circular 
binaries. The second time we allowed for initially circular orbits to see the effect of such 
systems on the circularization and cutoff periods. 

The distributions of circularization periods and cutoff periods resulting from the two 
Monte Carlo experiments described above allow us to determine the uncertainty on each of 
these measures. We will refer to the circularization period and the cutoff period derived 
from each binary population as the observed circularization period and the observed cutoff 
period. For reference two "parent" populations of ~ 20, 000 binary orbits were generated 
and tidally evolved. We will refer to the circularization period and the cutoff period derived 
from these parent populations as the true circularization period and the true cutoff period. 

Figure 5 show the distributions of the observed circularization periods (5a) and the 
observed cutoff periods (5b) derived for the 10,000 different binary populations excluding 
initially circular binary orbits. The true circularization period is shown in Figure 5a as 
a vertical solid line at 10.6 days. The distribution of circularization periods is slightly 
asymmetric with a tail at longer periods. Nonetheless, the true circularization period falls 
at the mode of the observed distribution. Thus, we take an observed circularization period 
as the best estimate of the true circularization period for any binary population. 

The interval enclosing 2/3 of the observed circularization periods on either side of the 
mode is marked by dotted vertical lines, and hereafter called the 2/3- interval. We define the 
uncertainty on the true circularization period of a binary population as the ± period intervals 
that the true circularization period can be shifted before the observed circularization period 
must be drawn from outside of the 2/3-interval. The errorbars on the true circularization 
period can thus be determined by considering two limiting cases: i) The observed circular- 
ization period is located at the short-period boundary of the 2/3 interval; 1.0 days shortward 
of the mode. In that case the observed circularization period would underestimate the true 
circularization period by 1.0 days. We assign a positive errorbar of 1.0 days to the true 
circularization period, ii) The observed circularization period is located at the long period 
boundary of the 2/3 interval; 1.2 days longward of the mode. In that case the observed 



^Orbits randomly selected with an initially negative eccentricity were excluded. 
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Fig. 5. — a) The distribution of circularization periods from 10,000 different binary pop- 
ulations each with 20 binary orbits with periods between 10 and 50 days. Only initially 
eccentric binary orbits (cj > 0.05) were considered. The vertical solid line mark the location 
of the true circularization period determined from a "parent" population of ~ 20,000 binary 
orbits (10.6 days). The interval containing 2/3 of all circularization periods around the mode 
is shown by the two vertical dotted lines at 9.6 and 11.8 days, b) The distribution of cutoff 
periods from the same 10,000 binary populations. The true cutoff period determined from 
the "parent" population of ~ 20,000 binary orbits is marked by a vertical solid line at 25.0 
days. The vertical dotted lines at 10.0 and 13.2 days mark the interval containing 2/3 of all 
cutoff periods around the mode at 11.3 days. 
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circularization period would overestimate the true circularization period by 1.2 days. We 
assign a negative errorbar of 1.2 days to the true circularization period. 

Accordingly, for a binary sample with ~ 20 orbits shorter than 50 days (e.g. the M35 
sample) the uncertainty range around the circularization period is + 1.0 days and - 1.2 days. 

The cutoff period of each binary population is defined here as the period of the longest 
period binary with an eccentricity below 0.05. The distribution of observed cutoff periods 
in Figure 5b is also slightly asymmetric with a tail toward longer periods and a mode of 
11.3 days. The true cutoff period is marked by a vertical sohd line at 25 days. Necessarily, 
the observed cutoff periods are shorter than the true cutoff period, and the mode of the 
distribution falls 13.7 days shortward thereof. 

The measurement uncertainty of the cutoff period is defined similarly to the circulariza- 
tion period. The 2/3- interval is marked by dotted vertical lines, and we define the uncertainty 
range around our estimate of the true cutoff period as the ± period intervals that the true 
cutoff period can be shifted before the observed cutoff period must be drawn from outside of 
the 2/3-interval. Accordingly, for a binary sample with 20 orbits shorter than 50 days (e.g. 
the M35 sample) the uncertainty range around the maximum likelihood estimate of the true 
cutoff period is -|- 1.3 days and - 1.9 days. 

We wish to emphasize that the maximum likelihood observed circularization period 
and the true circularization period arc the same, whereas the maximum likelihood observed 
cutoff period is offset from the true cutoff period by 13.7 days. The latter difference rcffects a 
significant sample size dependence of the cutoff period that obscures the physical dependence 
of the age of the binary population. 

In addition the Monte Carlo error analysis gives an error on the circularization period 
that is 1.0 day smaller than the error on the cutoff period. Thus we have shown that the 
circularization period is both more accurate and more precise than the cutoff period. 

Next we consider binary populations including initially circular (e < 0.05) binary orbits 
drawn from the Gaussian parent eccentricity distribution. Figure 6 shows the distributions 
of the observed circularization periods (6a) and the observed cutoff periods (6b) derived for 
10,000 binary populations. All vertical lines represent the same quantities as in Figure 5. 
Note the small effect of initially circular orbits on the distribution of the observed circular- 
ization periods and on the true circularization period derived from the parent population. 
The uncertainty range around the maximum likelihood estimate of the true circularization 
period is -|- 1.0 days and - 1.5 days, somewhat infiated by the initially circular orbits. 

In contrast, the distribution of cutoff periods is drastically changed by the presence of 
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Fig. 6. — a) The distribution of circularization periods from 10,000 different binary popu- 
lations each with 20 binary orbits with periods between 10 and 50 days. Initially circular 
binary orbits (cj < 0.05) were included. All vertical lines represent the same quantities as in 
Figure 5. The estimated true circularization period is 10.7 days, and the interval containing 
2/3 of all circularization periods around the mode range from 9.7 to 12.2 days, b) The 
distribution of cutoff periods from the same 10,000 binary populations. 
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initially circular binary orbits (e < 0.05). The true cutoff period is undefined and a tail of 
cutoff periods longward of the distribution mode is produced. 

In the presence of binaries with initially low eccentricity determination of the cutoff 
period is then dependent on the judgment of individual observers to define the eccentricity 
threshold between circular and eccentric systems and to select the binary defining the cutoff 
period. Examples of the effect of such judgments/selections on the value of the cutoff period 
will be given in Section 7 below. 

Comparison of Figure 5 and Figure 6 clearly demonstrates the robustness of the circular- 
ization period and the vulnerability of the cutoff period. Determination of the circularization 
period uses the information from all binary orbits, circular as well as eccentric, and is thus 
less vulnerable to initially circular orbits. The cutoff period is determined from one binary 
orbit alone and is thus vulnerable to the presence of initially circular orbits in the binary 
population, contamination from binaries with anomalous evolutionary paths, and/or binaries 
that have been falsely classified as cluster members. 

Importantly, the disadvantages of the cutoff period will be enhanced with larger binary 
populations, whereas the measurement uncertainty on the circularization period will decrease 
with increasing population size. 

To enable ourselves and colleagues to quote errors on circularization periods to binary 
populations of different sizes and ages, we repeated the Monte Carlo error analysis (allowing 
for initially circular binary orbits) 9 times. For each of 3 different sample sizes (10, 20, and 
30 binaries) we tidally evolved each of 10,000 samples so that the modes of the circularization 
period distributions were ~ 5, 10, and 15 days. Errors on the circularization periods were 
derived in the manner explained above. Table 2 list the derived errors and can be used as a 
reference for colleagues studying tidal circularization in clusters not listed in this paper. 

The coefficients (3 and 7 in the circularization function (eq. [1]) controls the steep- 
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ness/slope and the abruptness, respectively, of the function's transition from zero to non- 
zero eccentricities. Values of 0.14 for (3 and 1.0 for 7 were adopted to minimize the width 
of the distributions of circularization periods resulting from the Monte Carlo experiments 
described above, and to minimize the sensitivity of the circularization period to the choice 
of eccentricity threshold (e = 0.01) between circular and eccentric orbits. 

6. THE CIRCULARIZATION PERIOD IN M35 

After motivating, presenting and testing the circularization function and defining the 
circularization period, we are ready to apply this new diagnostic of tidal circularization to 
the population of binary orbits in M35. 

Figure 7 show the period-eccentricity distribution in the e — log(P) diagram of M35. The 
circularization function resulting in the minimum total absolute deviation between observed 
and calculated eccentricities is over-plotted and the locations of the observed circularization 
period (CP) and it's uncertainty interval are marked. The horizontal dotted line marks 
e = 0.01. The observed circularization period for M35 is 10.21};° days. The errors quoted 
were derived in Section 5 and arc listed in Table 2. 

The M35 period-eccentricity distribution show striking similarities to the artificially 
generated distributions created using eqs. (3) and (4) and an initial Gaussian eccetricity 
distribution. In addition to the period overlap between eccentric and circular orbits described 
in Section 3, 3 circular binaries with periods of ~ 10 days in M35 nicely represent the most 
frequent (typical) binaries located in the high density region at e < 0.05 and P ~ 10 days in 
Figure 4. The distribution of eccentricities of binaries with periods longer than ~ 10 days 
is consistent with a primordial Gaussian distribution and the circular binary (4037) with a 
period 16.49 days is likely to originate from the low eccentricity wing. 

Therefore, M35 is a good example of the advantages of the circularization period in 
defining the transition between eccentric and circular binary orbits. The best fit circulariza- 
tion function takes advantage of the information provided by all binary orbits and provide 
a more robust determination of the circularization period for the typical binary orbit. 

7. THE CIRCULARIZATION PERIODS OF 7 ADDITIONAL BINARY 

POPULATIONS 

Below we briefly discuss and present the period-eccentricity distributions and determine 
circularization periods for 7 additional binary populations spanning ages from ~ 3 Myr 
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Fig. 7. — The distribution of orbital eccentricity (e) as a function of log orbital period (P) for 
32 spectroscopic binary members of M35. The grey curve represent the best fit circularization 
function and the dotted horizontal fine mark e = 0.01. The circularization period (CP) is 
marked at log(P) = 1.0 (10.2 days) and it's uncertainty interval overplotted. 
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(PMS binaries) to ~ 10 Gyr (Galactic halo binaries). Figure 8a-h show the orbital data 
of each individual population with the best-fit circularization function over-plotted and the 
circularization period with error marked. All results are listed in Table 3; here we briefly 
discuss each population in turn. 

The PMS binary population (Figure 8a): The sample of 37 low-mass PMS binaries 
show the characteristic period overlap between eccentric and circular orbits. The orbital 
parameters are taken from Melo et al. (2001), and references to individual binaries can be 
found in their paper. The PMS sample is not strictly coeval, but cover an age range from ~ 1- 
10 Myr (Melo et al. 2001). We determine a circularization period of 7.11};^ days for the PMS 
binary sample. This value should be compared to the 7.56 day cutoff period determined by 
Melo et al. Note, that Melo et al. chose to disregard the circular orbit of binary RX J1301.0- 
7654a (P = 12.95 days (log(P) = 1.11), e = 0.04 ± 0.02; Covino, private communication) 
based on measurements of super-synchronous rotation of the individual binary components. 
Because the timescale for tidal synchronization is thought to be shorter than the timescale 
for circularization (Zahn 1977; Hut 1981), Melo et al. argue that the circular orbit of RX 
J1301.0-7654a might not be a result of tidal circularization. However, as mentioned by Melo 
et al., super-synchronous stellar rotation is an expected result of the transition from the PMS 
to the Z AMS due to conservation of angular momentum and a decline in the efficiency of tidal 
breaking as the stars contract. Thus we have chosen to include binary RX J1301.0- 7654a. 
However, due to the robust nature of our new diagnostic, excluding binary RX J1301.0- 7654a 
from the PMS sample causes only a 0.5 day decrease in the circularization period, whereas 
the longest period circular orbit would decrease by 5.39 days (from 12.95 to 7.56 days). 

The Pleiades binary population (Figure 8b): The eccentricity distribution of 12 binary 
orbits in the ~ 100 Myr Pleiades cluster provides poor samphng of the transition region. Also, 
unlike any other population included in this study, the Pleiades features a gap between the 
longest period circular and the shortest period eccentric orbit. By fitting the circularization 
function we make use of the information provided by the binary orbits on both sides of the 
gap, and determine a circularization period of 7.2l}'g days. This value should be compared 
to the 7.7 day cutoff period previously used for the Pleiades cluster. The orbital parameters 
shown in Figure 8b are taken from Mermilliod et al. (1997, 1992). The estimated masses for 
the binary primary components falls within ~ 0.9 Mq — 1.4 Mq. 

The M35 binary population (Figure 8c): The distribution of orbital eccentricities in 
M35 (discussed above) is shown here for the purpose of comparison. 

The Hyades/Praesepe binary populations (Figure 8d): The sample of 47 binary orbits 
in the ~ 630 Myr Hyades/Praesepe twin clusters show 3 eccentric binaries shortward of the 
longest period circular orbit. The best fit circularization function is over-plotted and the 
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Fig. 8. — The period-eccentricity distributions for 8 late-type binary populations: a) PMS, 
b) Pleiades, c) M35, d) Hyades/Praesepe, e) M67, f) NGC188, g) field, and h) halo. Over- 
plotted each distribution is the best fit circularization function. A solid horizontal line mark 
e(P) = 0.00. The circularization period (CP) and it's uncertainty are marked on the period 
axis by a vertical black line and a grey bar, respectively. The circularization period and the 
age of the binary population are printed in each plot. 
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circularization period of 3.2l]^ J days is marked. The cutoff period for the Hyades/Praesepe 
sample has previously been set by the 8.49 day circular binary J331. J331 consist of two 0.5 
Mq components, and thus is substantially different from the binaries with approximately 
solar-mass components defining the transition regions in other binary populations. For 
this reason Mathieu et al. (1992) questioned whether binary J331 should be disregarded, 
as the difference in mass and stellar structure (depth of convective envelope) might have 
significantly altered the timescale of tidal circularization. We will not comment further on 
the effect of stellar mass here, but simply emphasize that while we have included binary J331, 
disregarding it leads to no noticeable change in the Hyades/Praesepe circularization period, 
compared to a -2.75 days change from 8.49 to 5.74 days in the Hyades/Praesepe cutoff period. 
The orbital parameters shown in Figure 8d are taken from Griffin et al. (1985, 1982); Griffin 
& Gunn (1981, 1978) for Hyadcs binaries and Mermilliod & Mayor (1999); Mermilliod ct al. 
(1992, 1990) for Praesepe binaries. The estimated masses for the the primary components 
of this sample falls within ~ 0.5 Mq — 1.5 Mq. 

The M67 binary populations (Figure 8e): The sample of spectroscopic binaries in the 
old open clusters M67 contain stellar components ranging in evolution from the unevolved 
main-sequence to the tip of the giant branch (Mathieu, priv. comm.). As the rate of tidal 
circularization depends sensitively on stellar radius and the depth of the convection zone, 
we will consider only binaries in M67 with unevolved primary components. Figure 8e shows 
the distribution of orbital eccentricity vs. log orbital period for 39 unevolved binaries in 
M67. Two eccentric orbits are found shortward of the longest period circular orbit. The 
best fit circularization function is over-plotted and the circularization period of 12.llJ;5 days 
is marked. The mass-range of the binary primary components is ~ 0.9 Mq — 1.2 Mq. 

The NGC188 binary populations (Figure 8f): As for M67, the binary population in 
NGC188 contains evolved as well as unevolved stellar components. Again we will consider 
only binaries with unevolved primary components. The distribution of orbital eccentricity 
vs. log orbital period for 27 unevolved binary stars in NGC188 is shown in Figure 8f. The 
binary population contains a single eccentric system with a period shorter than the longest 
period circular binary. The circularization period of 14.5^2 2 days is marked and the best 
fit circularization function is over-plotted. The orbital data shown in Figure 8f are taken 
from Mathieu et al. (2004). The mass-range of the primary components in this sample is 

~ 0.95 Mq - 1.05 Mq. 

The field binary population (Figure 8g): The orbital characteristics of 50 close solar-type 
binaries from the solar neighborhood (hereinafter: the field) is shown with the circularization 
period of 10.3^3 '^ days marked and the best fit circularization function over-plotted. The 
orbital data are taken from Duquennoy &: Mayor (1991), who adopt an age range of 7-11 Gyr 
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for the population of field binaries. We adopt 9 Gyr as the age of the field sample. However, 
the age of some individual binary systems are poorly known, thus despite the well-defined 
circularization period, the age-ambiguities of this sample makes it less reliable for probing 
the evolution of tidal circularization. 

The Galactic halo binary population (Figure 8h): The orbital characteristics of 61 
binaries from the Galactic halo. The best fit circularization function is over-plotted and the 
circularization period of 15.6^3 2 days is marked. The orbital data are taken from Latham 
et al. (2002). We adopt an age of 10 Gyr for the halo binary population. From a large 
sample of 171 high proper motion spectroscopic binaries, Latham ct al. (2002) find no obvious 
differences between the binary characteristics in the halo and in the disk populations. The 
observed frequency is the same and the period distributions are consistent with the hypothesis 
that the binaries are drawn from the same parent population. Nonetheless, the halo binary 
population is not strictly coeval and the stars differ significantly in mass and metallicity 
from those of the open cluster/disk populations. Arguably, comparison of solar-mass/solar- 
met allicity tidal circularization theory to the halo cutoff period requires consideration of 
these differences. 

It is important to note that in the majority of the binary populations discussed here, cir- 
cular or marginally circular binary orbits are observed with periods from tens to hundreds of 
days longer than the circularization period. The selection of unevolved main-sequence bina- 
ries for each population suggest that these long period (P ~ 10"*^ — 10^ days) low eccentricity 
(e ~ 0.05) orbits are not a result of accelerated tidal circularization due to evolved primary 
components. The existence of long-period low eccentricity binaries in the PMS, Pleiades and 
M35 samples is a strong indication that binary stars can form with circular or marginally 
circular orbits, consistent with an initial Gaussian distribution of orbital eccentricities. 

The binary populations shown in Figure 8 also show a high frequency of high-eccentricity 
orbits at periods shorter than the one or more circular orbits. The ability to reproduce such 
systems from binaries in the high-eccentricity tail of the initial Gaussian distribution by 
numerically integrating current theoretical models of tidal circularization (Duquennoy et al. 
1992, Section 5 this paper), provide further support for the existence of an initial Gaussian 
distribution of orbital eccentricities. 



8. EVOLUTION OF TIDAL CIRCULARIZATION 

While different tidal circularization theories agree that the timescale of circularization 
depends strongly on stellar separation, and that consequently a transition from eccentric 
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to circular orbits is expected at a well-defined binary period, the different theories describe 
different mechanisms for the tidal dissipation leading to circularization. The mechanism and 
efficiency of tidal dissipation in a given model reflects itself in the model's prediction of the 
rate of circularization and the evolution of the tidal circularization period. Observational 
determinations of the circularization periods in binary populations of different ages provide 
a critical test of the dissipation mechanism(s) responsible. 

Detailed discussion of the distribution of circularization cutoff periods with age, and its 
ability to constrain tidal circularization theory, was given by Witte & Savonije (2002), Melo 
et al. (2001), and Mathicu et al. (1992). However, the addition of the new large binary 
populations of M35 and NGC188 (Mathieu et al. 2004), the introduction of a new robust 
measure of tidal circularization (the circularization period) with measured uncertainties, and 
the recent theoretical work on dynamical tides in late- type stars (e.g. Witte & Savonije 2002), 
motivate a fresh look at how the theory compares to the observations. 

We show in Figure 9 the circularization periods for all binary populations discussed in 
this paper (solid/open circles) at their respective ages. Errorbars show the uncertainties on 
the circularization periods based on the Monte Carlo experiments described in Section 5. 

With the exception of the Hyades/Praesepe population, the circularization periods of 
the coeval solar-type binary populations show a steady increase with age from the PMS and 
early main-sequence to the late main-sequence phase. The circularization period of the halo 
and field binaries follow this trend. However, due to the non-coevality and uncertainty in 
age of the field sample, and the differences in mass and metaUicity of the stars in the halo 
sample, we choose not to include their circularization periods in the subsequent analysis. 

The small value of the circularization period of the Hyades/Praesepe sample can be 
explained in part by the high number of short period eccentric binaries, in particular the 
two highly eccentric short period binaries (Hyadcs: HD30738, P = 5.75 days, e — 0.354; 
Praesepe: KW181, P = 5.87 days, e = 0.357). However, even in the absence of these two 
systems the Hyades/Praesepe circularization period will be 5.9 days and thus still deviate 
from the overall trend of increasing circularization periods with increasing age. We note that 
Duquennoy & Mayor (1991) draw attention to an excess of short period binaries (log(P) < 
1 day) in the Hyades period distribution when compared to distribution from their G-dwarf 
sample. 

The predictions of main-sequence tidal circularization in the framework of both the 
equilibrium tide theory and the dynamical tide theory are also shown in Figure 9, as is the 
prediction by Zahn & Bouchet (1989) including PMS tidal circularization. Note that these 
different theoretical predictions displayed in Figure 9 use different criteria for determination 
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Fig. 9. — The distribution of circularization periods with age for 8 late-type binary popula- 
tion (solid/open circles). Errorbars represent the uncertainties on the circularization periods 
derived in Section 5. The solid curve shows the predicted cutoff period as a function of time 
based on main-sequence tidal circularization using the revised equilibrium tide theory by 
Zahn (1989) (Claret & Cunha 1997). The broad dashed band represent the predicted cutoff 
period period for initially super-synchronous 1 Mq stars calculated in the framework of the 
dynamical tide model including resonance locking (Witte & Savonije 2002). The horizontal 
grey band represents the prediction by Zahn & Bouchet (1989) in which tidal circularization 
is significant only during the PMS phase. 
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of the cutoff period. Claret & Cunha (1997, equilibrium tide theory) define the cutoff 
period at a given age as the longest orbital period for which the relative variation in the 
eccentricity is 0.5% of the initial value (ej„j = 0.30). Witte & Savonije (2002, dynamical tide 
theory) define the cutoff period as the longest orbital period for which a binary with 1 Mq 
components has been circularized to e < 0.01. As such the predictions of Witte & Sovonije 
are upper limits on the circularization period. Zahn & Bouchet (1989, PMS tidal theory) 
define the cutoff period as the period at which a binary with component masses between 
0.5 Mq — 1.25 Mq and initial orbital eccentricity of 0.2 or 0.3 circularizes to e < 0.005. 
Recall that the circularization periods displayed in Figure 9 represent the orbital period at 
which a binary orbit with initial eccentricities in the range 0.2 - 0.5 evolves to e = 0.01 at 
the age of the population. 

The equilibrium tide theory: Let us consider first the equilibrium tide theory using the 
dissipation mechanism by Zahn (1989). The equilibrium tide theory predicts active tidal 
circularization throughout the main-sequence phase, leading to longer cutoff periods with 
increasing population age (solid curve in Figure 9). However, as pointed out by Claret & 
Cunha (1997), an artificial enhancement of the turbulent dissipation is necessary to fit the 
observed cutoff periods. Comparison with the new circularization periods, again with the 
exception of the circularization period of the Hyades/Praesepe population, show a similar 
discrepancy between the observations and the theoretical predictions of pure main-sequence 
tidal circularization. 

The grey horizontal band in Figure 9 shows the prediction by Zahn & Bouchet (1989) of 
tidal circularization including the PMS phase. Based on the reduced rate of tidal circulariza- 
tion derived by Zahn (1989), Zahn & Bouchet suggested that all tidal circularization occurs 
during the PMS phase when the stars have larger radii and deeper convective envelopes. 
They found that for main-sequence stars the efficiency of turbulent dissipation derived from 
the equilibrium tide theory is so small that tidal circularization following the PMS phase 
is negligible despite the Gyr timescales available. Their PMS tidal circularization theory 
predicts a range of cutoff periods between ~ 7.2 - 8.5 days for binaries with components 
with masses between 0.5 - 1.25 Mq and initial eccentricities of either 0.2 or 0.3. Accordingly, 
all late-type main-sequence binaries should have circular orbits for period less than ~ 8 days 
and (primordial) eccentric orbits at longer periods, independent of age. 

The hypothesis that PMS tidal circularization alone sets a tidal cutoff period, inde- 
pendent of age, is not consistent with the observed distribution of circularization periods. 
We performed a t^st fitting different constant theoretical circularization periods to the 
distribution of observed circularization periods. We find that there is less than a 1% prob- 
ability that the observed circularization periods (excluding or including the field and halo 
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samples) derive from a model predicting no evolution of the circularization periods with age. 
Repeating the test for all binary populations but the Hyades/Praesepc population also 
tells us that there is less than 1% chance that the PMS, Pleiades, M35, M67, and NGC188 
populations derive from a model predicting pure PMS circularization. 

The dynamical tide theory. The predicted evolution of main-sequence tidal circulariza- 
tion using the dynamical tide theory with inclusion of resonance locking (Wittc & Savonije 
2002) is shown in Figure 9 as a broad dashed band. The width of the dashed band represents 
the range of predicted evolutions of the cutoff period for binaries with stellar components 
rotating at super- synchronous {Q = 2ujp), pseudo-synchronous {Q = ujp), and orbital (f2 = uj) 
speeds. Q, LUp, and u denote the initial stellar angular velocity, the orbital angular velocity at 
periastron, and the average orbital angular velocity, respectively. The masses of the binary 
components used by Witte & Savonije are 1 Mq. 

Regardless of stellar rotation velocity, the dynamical tide theory with inclusion of the 
resonance locking mechanism (Witte & Savonije 2002) appears to predicts shghtly more 
efficient tidal circularization than the equilibrium tide theory throughout the main-sequence 
phase. However, Witte & Savonije (2002) define the cutoff period as the longest orbital 
period for which a binary has been circularized to e < 0.01. This definition implies that the 
predicted cutoff periods are likely set by binaries with low initial orbital eccentricities, and 
thus are not directly comparable to the predictions by Claret & Cunha (1997). Still there is 
a discrepancy between the predicted level of circularization and the distribution of observed 
circularization periods. For the special case of very slowly rotating stars, the calculations 
of Witte & Savonije (2002, not shown in Figure 9) offer a near match to the circularization 
periods of Pleiades, M35, M67 and NGC188. However, the required stellar rotation periods 
are of order 100 days for stars in close binary systems with ages between 0.1 and 10 Gyr. Such 
very slow rotations would be unexpected and thus needs to be observationally established. 

The slopes of the main-sequence tidal circularization models match the increase in cir- 
cularization period seen in Figure 9. Thus if the efficiency of tidal dissipation were to be 
artificially enhancement for both the equilibrium tide and the dynamical tide theories, the 
predicted main-sequence evolution for both theories would be in agreement with the observed 
increase in circularization period for all populations but the Hyades/Praesepe. 

The hybrid scenario: The lack of success of both PMS and main-sequence tidal cir- 
cularization theories to account for the observations of tidal circularization motivated the 
suggestion of a "hybrid scenario" (Mathieu et al. 1992). The hybrid scenario is heuristic 
and was motivated by a distribution of cutoff periods that appeared age-independent for 
populations younger than ~ 1 Gyr while showing a positive correlation with age for popu- 
lations older than ~ 1 Gyr. The hybrid scenario suggest that tidal circularization in binary 
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populations younger than ~ 1 Gyr derive from PMS tidal circularization, and that after the 
passage of ~ 1 Gyr the integrated main-sequence tidal circularization begins to circularize 
binaries with orbital periods of ~ 10 days. 

Considering the new distribution of circularization periods the need for a hybrid scenario 
is uncertain. With the exception of the Hyades/Praesepe population the observed circular- 
ization periods increase with age from the PMS and throughout the main-sequence phase. 
A test shows that there is less than 1% probability that the circularization periods for 
populations younger than ~ 1 Gyr derive from a model predicting no main-sequence tidal 
circularization. Excluding the Hyades/Pracscpc sample result in value of ~ 1. Thus, 
the current distribution of circularization periods for the youngest binary populations can- 
not distinguish between a model predicting an age-independent circularization period or a 
model predicting a continuous increase in the circularization period from the PMS onward. 

If tidal circularization is in fact active throughout the main-sequence phase, the ques- 
tion remains which if any of the suggested dissipation mechanisms is responsible. Either the 
equilibrium or the dynamical tide mechanisms, or perhaps both in combination, could be 
responsible for the circularization of binary orbits during the main-sequence phase. How- 
ever, unless both mechanisms are at work simultaneously, more efficient tidal dissipation is 
needed in both the equilibrium and the dynamical tide theories. We note that PMS tidal 
circularization is still needed to explain the 7.1 day circularization period of the PMS sample, 
and that an increase in the efficiency of turbulent dissipation of the equilibrium tide theory 
presumably will affect circularization during the PMS as well. 

Alternatively, assuming pure PMS tidal circularization, the differences in circularization 
periods among the binary populations in Figure 9 could be due to differences in initial 
stellar and circumstellar conditions. Goodman & Dickson (1998) speculate that perhaps the 
differences in mass and/or metallicity between the younger disk binaries and the older disk 
and halo binaries could be the cause of differences in cutoff periods set during the PMS phase. 
While the ranges of primary stellar masses and metallicities of the populations included in 
this study are similar, other stellar/circumstellar and binary /circumbinary parameters might 
be of importance to the efficiency and duration of tidal circularization during the PMS phase. 
The strong dependence of the efficiency of circularization on stellar rotation displayed in the 
dynamical tide theory by Witte & Savonije (2002) is a powerful example of how differences 
in initial conditions may affect the evolution of tidal circularization. 

It is difficult to draw secure conclusions about the evolution of tidal circularization 
based on only 6 circularization periods. The future observational goal is to further populate 
the age vs. circularization period diagram (Figure 9) with reliable circularization periods 
at carefully selected ages based on large samples of binary orbits from homogeneous coeval 
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binary populations. Following the present study and the recent publication of binary orbits in 
NGC188 (Mathieu et al. 2004), we plan to contribute (work in progress) such circularization 
periods at the ages of the young open cluster M34 (~ 250 Myr) and the intermediate age 
open cluster NGC6819 (~ 2.5 Gyr). By inspection of Figure 9 it is evident that these future 
circularization periods will greatly improve the observational constraint on future theoretical 
models and predictions of tidal circularization in close late-type binary stars. 

The future goal for theory of tidal circularization in late-type main-sequence stars is 
to construct models that combine PMS tidal circularization with main-sequence tidal cir- 
cularization, and predict more efficient tidal circularization in binaries with main-sequence 
components to account for the observed evolution of circularization. 

9. SUMMARY AND CONCLUSIONS 

In a coeval population of binary stars the closest binaries tend to have circular orbits, 
while wider binaries have orbits with a distribution of non-zero eccentricities. Determination 
of the orbital period to which binaries have been circularized is an important observational 
constraint on theoretical models seeking to explain the dissipation mechanism(s) responsible 
for tidal circularization. Coeval populations of short-period binary orbits are therefore key 
observational contributions. 

However, despite homogeneity in age, mass and metallicity, a population of short-period 
binaries may be small in size and will consist of systems with a distribution of initial ec- 
centricities. The size of a binary population will affect the uncertainty on its circularization 
period. The initial eccentricity of a circular binary with a period equal to the circularization 
period is unknown, but can be estimated through the use of theoretical models. Proper 
comparison between observed and theoretical circularization periods thus require measure- 
ment uncertainties on the observed circularization periods and information about the initial 
eccentricities of the binaries that circularize at those periods. 

The primary contributions of this study are: 1) a new population of solar-type spec- 
troscopic binary orbits for the 150 Myr open cluster M35; 2) a new diagnostic to determine 
the period of circularization of the most frequently occurring binary; 3) an evaluation of 
the performance of this new diagnostic providing measurement uncertainties on the result- 
ing circularization periods and information about the initial eccentricity of the binaries that 
circularize at that period; and 4) a comparison of the distribution of circularization periods 
from 8 binary populations to the predictions from theoretical models. 

The new sample of 32 spectroscopic binary orbits in the "Pleiades age" (~ 150 Myr) 
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open cluster M35 greatly improves the constraint on tidal circularization at the early stage 
of the main-sequence phase. Prior to this study, the small sample of binary orbits in the 
Pleiades provided only very rough limits on the transition between eccentric and circular 
orbits at this age. 

Monte Carlo simulations of tidal circularization using the equilibrium tide theory by 
Zahn (1977) were used to create artificial e— log(P) distributions. These simulations, together 
with the e — log(P) distributions of M35 and 7 additional published binary populations, were 
used to critically assess the adequacy of the tidal circularization cutoff period (the period 
of the longest period circular binary) to measure the degree of tidal circularization. We 
conclude that the cutoff period is not an optimal measure of the transition between eccentric 
and circular orbits because of the following disadvantages: 

• The simulations of tidal circularization using the equilibrium tide theory by Zahn 
(1977) show that the cutoff period is hkely to originate from the low eccentricity tail 
of an initial Gaussian eccentricity distribution. Therefore, the cutoff period does not 
represent the period of circularization of the most frequently occurring binary orbit. 

• The expected value of the observed cutoff period of a population of binaries varies with 
the size of the population. 

• The cutoff period is vulnerable to the presence of circular orbits that are not a measure 
of tidal effects, such as primordial circular binary orbits or orbits of binaries with 
anomalous evolutionary paths. 

• Determination of the cutoff period makes no use of the information provided by eccen- 
tric orbits in the e — log(P) diagram. 

Motivated by the need for a more robust determination of the orbital period at which 
a binary orbit with the most frequent initial orbital eccentricity circularizes, we introduce 
the circularization function and the circularization period. The functional shape of the 
circularization function (eq. [1]) is fixed and mimics the transition from eccentric to circular 
orbits in the e — log(P) diagram (Figure 4). The circularization period is defined as the 
period for which the circularization function has a value of 0.01 (e(CP) = 0.01). The 
location in period of the circularization function thus determines the circularization period. 
The advantages of the circularization function/period are: 

• Numerical simulations of tidal circularization using the equilibrium tide theory by Zahn 
(1977) show that the circularization period is the period at which a binary with the 
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most frequent initial eccentricity (e ~ 0.35) evolves to e = 0.01 at the age of the 
population. 

• The circularization period is determined from all binary orbits (circular and eccentric) , 
and so is less vulnerable to primordial circular binaries or binaries with anomalous 
evolutionary paths. 

• The observed circularization period is the best estimate of the true circularization 
period independent of the size of the binary population. 

Monte Carlo error analysis allows us to examine the distribution of the circularization 
periods. Prom this analysis we estimate and provide the uncertainty on the circularization 
period as a function of the size of the binary sample and of the value of the circularization 
period. We note that even in the absence of initial circular binaries and/or contamination 
from binaries with anomalous evolutionary paths, the Monte Carlo error analysis show that 
the uncertainty on the circularization period is smaller than the uncertainty on the cutoff 
period. Thus we conclude that the circularization period is both more accurate and more 
precise than the cutoff period. 

We present the circularization periods for M35 and 7 additional binary populations 
ranging in age from ~ 3 Myr to ~ 10 Gyr. With the exception of the circularization period 
of the Hyades/Praesepe population, the distribution of circularization periods with age show 
a steady increase from the PMS to the late main-sequence. 

We compare the distribution of circularization periods to the predictions by theoretical 
models. The models of main-sequence tidal circularization using either the equilibrium tide 
theory (Zahn 1989) or the dynamical tide theory with resonance locking (Witte & Savonije 
2002) both predict longer cutoff periods with increasing population age in agreement with 
the trend in the distribution of circularization periods. However, the predicted circular- 
ization periods fall significantly below the observed circuleirization periods. This 
suggest that the efficiency of the dissipation in these models is too low. 

A model including PMS tidal circularization (Zahn & Bouchet 1989) predicts negligible 
circularization during the main-sequence phase. This prediction is not consistent with the 
observed distribution of circularization periods. 

The need for a hybrid scenario as proposed by (Mathieu et al. 1992) is uncertain. The 
current distribution of circularization periods for the youngest binary populations cannot 
distinguish between a model predicting an age-independent circularization period or a model 
predicting a continuous increase in the circularization period from the PMS onward. 
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The goal for future theories of tidal circularization in late-type main-sequence binaries 
is to combine PMS and main-sequence circularization in one model that can account for the 
circularization periods observed at all ages. 
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Table 3. Distribution of Circularization Periods with Population Age 
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